An eccentric disk model for the nucleus of M31 
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ABSTRACT 

The nucleus of MSI may be a thick eccentric disk, composed of stars 
traveling on nearly Keplerian orbits around a black hole or other dark 
compact object. This hypothesis reproduces most of the features seen 
in HST photometry of the center of MSI; in particular the bright off- 
center source PI is the apoapsis of the disk. An eccentric disk can also 
explain the rotation curve and asymmetric dispersion profile revealed 
by ground-based observations. The central object must be smaller 
than ~ 1 pc so that the potential felt by the disk is nearly Keplerian. 
The disk eccentricity may be excited by a secular instability driven 
by dynamical friction from the bulge. 
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1 Introduction 


The nucleus of M31 may be the nearest dead quasar, that is, it may contain a black 
hole (hereafter BH) of sufficient mass (~ 10® Mq) to power an active galactic nucleus 
when offered an adequate fuel supply. The nucleus was first resolved in high-resolution 
(FWHM=0”2) photographs by the Stratoscope II balloon-borne telescope (Light, Daniel¬ 
son & Schwarzschild 1974), which also revealed that the peak brightness of the nucleus 
was offset relative to its outer parts. 

Hubble Space Telescope (HST) images with FWHM=0”08 show that the nucleus 
contains two separate components, separated by 0”49 (Lauer et al. 1993, hereafter L93; 
see also King, Stanford & Crane 1995). The component with the lower surface brightness, 
P2, coincides with the center of the outer isophotes of the nucleus and the bulge photocenter 
to within ~ 0”05, and exhibits a shallow central cusp. The brighter component, PI, is 
well-resolved and coincides with the nuclear core measured by Stratoscope. No strong color 
gradients are visible between the near-infrared and the far ultraviolet, except for enhanced 
emission in the far ultraviolet close to P2 (Mould et al. 1989; King et al. 1995). 

Ground-based spectroscopy (Dressier & Richstone 1988; Kormendy 1988; Bacon et al. 
1994; van der Marel et al. 1994) shows that the rotation speed rises from near zero outside 
4” to a maximum of ~ 120kms“^ at 1 " from the center. The location and amplitude of 
the rotation maximum are determined by the spatial resolution, so the actual maximum 
is almost certainly larger. The center of symmetry of the rotation curve coincides with 
P2 to within 0”2; the curve is symmetric so there is no evidence that the brighter peak 
PI perturbs the rotation curve. The velocity dispersion rises from ~ 150kms“^ outside 
2 " to a maximum of 220kms“^; once again the actual maximum is certainly larger. The 
dispersion peak is displaced from P2 to the anti-Pl side and the dispersion profile is not 
symmetric. 

Kinematic models—which usually neglect these asymmetries—require that the mass 
in the central arcsecond is dominated by a BH or other compact non-luminous mass con¬ 
centration, with mass M, ~ 10^'^ Mq (Kormendy 1988; Dressier & Richstone 1988; Bacon 
et al. 1994). 

Explanations of the asymmetries in the M31 nucleus are discussed by L93 and Bacon 
et al. (1994). Partial obscuration by dust is implausible because there are no color gradients 
or re-radiated far-infrared emission from the dust; also, if P2 is obscured by an asymmetric 
dust distribution, why should it be exactly at the center of the bulge? Bacon et al. (1994) 
suggested that PI is off-center because of a dynamical oscillation, but no specific model 
was proposed so their suggestion cannot be assessed. 

A natural hypothesis is that PI is a separate stellar system orbiting the nucleus. 
However, assuming PI is close to the nucleus—the probability of a chance superposition 
of a foreground or background object is negligible—its orbit rapidly decays by dynamical 
friction from the bulge. For the best-fit bulge model described below, the frictional decay 
time^ for an object of mass m is ~ 6 X 10^yr(10® Mq/ht), for radii between 0?5 and 3” 

^ The object is assumed to be on a circular orbit and the decay time is L/L where L is the angular 
momentum. The velocity distribution of bulge stars is assumed to be Maxwellian and the drag is computed 
using the usual Chandrasekhar formula with Coulomb logarithm In A = ln(3.0). 
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from the center (for comparison the orbital time at 1" is 8 X 10"^ yr). Thus the present 
configuration is short-lived. Most attempts to evade this conclusion are implausible: de¬ 
pletion of the stars that contribute to the friction does not occur because tidal forces in the 
triaxial bulge (Stark 1977) continually bring fresh stars close to the center; compensating 
the friction by energy transfer from a disk is only a temporary fix, because the reservoir 
of energy in the disk is too small. Finally, tidal forces from the central BH should disrupt 
PI unless its mass-to-light ratio is substantially higher than a normal stellar population. 

Some of these problems are avoided if PI is a massive black hole surrounded by bound 
stars. Lacey & Ostriker (1985) have argued that the dark matter in galactic halos may 
be composed of BHs of mass ~ 10®'^ Mq (but see Moore 1993 and Carr 1994 for counter¬ 
arguments). Black holes of this mass will spiral to the center by dynamical friction from 
as far out as ~ 1 kpc, with one arriving every 10® yr or so. The fate of these holes is 
controversial: (a) they may all coalesce into a dense cluster or even a single BH (Hut & 
Rees 1992); (b) gravitational interactions between a binary BH and a single BH may eject 
one or both systems, so that the number of BHs at the center is normally 0, 1, or 2 but 
no higher (Xu & Ostriker 1994). Thus the short lifetime of PI is not a problem, because 
fresh black holes arrive at a steady rate; in fact, Xu and Ostriker show that a binary BH 
with separation similar to that of PI and P2 is present about 10% of the time. However, 
it remains uncertain whether this model can naturally explain the coincidence of P2 with 
the center of the bulge and the symmetry of the rotation curve; also, acquiring the “cloak” 
of stars around the BH at PI is difficult, since the phase-space density of stars around PI 
is higher than anywhere else in the bulge or nucleus. 

The purpose of this paper is to present a model for the central region of M31 that 
explains the observed asymmetries. 


2 A model for the nucleus of M31 

The model assumes that a single massive BH is located at the center of the bulge, which 
may be taken to be the center of P2 (apart from a small correction described below). 

The nucleus is assumed to be an eccentric stellar disk, composed of stars traveling on 
Keplerian orbits around the BH at P2. The disk is eccentric because the apsides of the 
orbits are approximately aligned. Kormendy (1988) has already argued that the nucleus 
may be a disk; the only novelty here is that the disk is eccentric. 

The goal of this section is to demonstrate that an eccentric disk can explain the main 
features of the M31 nucleus; the model is not designed to fit the observations to the highest 
possible accuracy. The distance to M31 is assumed to be 770 kpc (Freedman & Madore 
1990) so that 1" = 3.73 pc, and the absorption correction is taken to be Ay = 0T24. 
The mass of the BH is M. = 7.5 X 10^ Mq and the mass-to-light ratio of the bulge is 
M/L = 4 M 0 /L 0 after correcting for absorption, values chosen to yield approximate 
agreement with the observed dispersion profile (see below). The nuclear stars are assumed 
to have the same mass-to-light ratio as the bulge. The mass-to-light ratio of the bulge is 
close to the estimate 5.6 M 0 /L 0 given by Kent (1989, after correcting for absorption and 
converting from r to V). 
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The bulge is assumed to be spherical (its actual ellipticity e ~ 0.15-0.20; Kent 1989) 
and its surface brightness is described by the fitting formula 
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which reproduces the surface brightness distribution of a wide variety of elliptical galaxies 
and bulges (Lauer et ah 1995). The formula has five parameters: the three exponents a, /3, 
7 , the break radius and the surface brightness at the break radius, Rotation of the 
bulge is neglected. The velocity dispersion tensor in the bulge is assumed to be isotropic, 
and is determined by solving the Jeans equations in the combined potential field of the 
central BH and the bulge stars. 

The nucleus (or “nuclear disk”) is an eccentric disk lying in the plane defined by the 
disk of M31 (inclination 77°, position angle of the line of nodes 38°). The disk is represented 
by assuming that its stars he on a small number of nested Keplerian orbits (“ringlets”) 
with semi-major axes a„, eccentricities e„, and total luminosities n = 1,. . . ,iV. The 
focus of each ringlet is at P2 and the apoapsis direction is assumed to coincide with the 
P2-P1 axis (position angle 43° according to L93, which implies that the azimuthal angle 
in the disk plane between the line of apsides and the line of nodes is 21°). To smooth 
the model and to account crudely for the disk thickness, the sky brightness distribution of 
each ringlet is convolved with a Gaussian point-spread function with standard deviation 
cr„; thus cr„/a„ is an approximate measure of the eccentricity and inclination spread in the 
disk. The velocity field is determined from the Keplerian velocities of the ringlets, with no 
correction for asymmetric drift. 

The model neglects the gravitational influence of the nuclear disk on the bulge, since 
its mass is only 16% of the BH mass in the best-fit modeP; the model also neglects the 
influence of the bulge on the nuclear disk, since its mass interior to 1" is only 2% of the 
BH mass. 

The parameters of the model are determined by fitting two data sets: 

(a) The surface brightness measured by HST in a strip of width 0”22 along the P2-P1 
axis, extending to ±4” from P2 (see Figure 6 of L93). 

(b) The surface brightness on the major axis of the bulge between A" and 100” (Kent 
1989). Kent’s r-band magnitudes are increased by 0T36 so that they match the HST 
K-band magnitudes in the region of overlap. 

The fitting procedure minimizes the mean square deviation in magnitude, giving equal 
weight to the two data sets. The model shown in this paper (Figure 1) uses only 3 ringlets, 
but already fits the data reasonably well (the rms residual in either data set is less than 
0T04); models with more ringlets provided only modest improvements to the fit and are 
increasingly ill-conditioned. The model does not reproduce precisely the shape of the profile 
in the region between P2 at x = 0 and PI near x = —0”6; however, these discrepancies 
could probably be removed in more careful models and are nowhere greater than OTl. 


® The principal consequence of this approximation is that our estimate of the BH mass is probably too 
large by ~ 10%. 
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1. (a) V-band surface brightness in a slice 0!^22 wide along the P2-P1 axis. P2 is at x = 0 and 

PI is at X = —0!^5. The hlled circles show the HST data (Figure 6 of L93), and the solid line 

shows the best-ht model with 3 ringlets. The dashed and dotted lines show the contributions 

from the bulge and the nuclear disk, (b) V-band surface brightness along the major axis of the 
bulge. The hlled circles show r-band photometry from Kent (1989), adjusted by 0T36 downward 
to agree with L93 in the region of overlap. The two solid lines show the best-ht model along both 
sides of the P2-P1 axis; the dashed and dotted lines show the contributions from the bulge and 
the nuclear disk. 

The fit shown in Figure 1 corresponds to the parameters (luminosities are corrected 
for absorption but apparent magnitudes are not), 

a = 0.58, /3 = 1.37, 7 = 0.15, /& = 17”42 (sec)“^ n = A9'.'l, 

ai = 0"49, ei = 0.44, Li = 5.8 x 10® L© = 15”09, ai = 0"16, 

02 = 0"87, 62 = 0.15, L 2 = 9.5 X 10® L© = 14”56, 02 = 0"29, (2) 

03 = 1"53, 63 = 0.04, Ls = 1.43 x 10® L© = 14”12, 03 = 0"78. 

The total luminosity of the disk is 3.0 X 10® Lq = 13T3, although the peak at PI is mostly 
due to ringlet 1 which is fainter. Many of these parameters are rather uncertain as the 
inversion is ill-conditioned. 

Figure 2 shows the surface brightness distribution predicted by the model, which can 
be compared to Figure 2 of L93. The model successfully reproduces the main features of 
the HST image, even though the fit was based on data from only the P2-P1 axis. There 
are minor discrepancies—for example, the ellipticity of the isophotes surrounding PI and 
P2 is higher in the model than in the data—but these should be easy to remove in more 
careful models. 
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2 . Contour map of the surface brightness of the best-ht model. The contour interval is 0™1. 
The map can be compared to Figure 2 of L93. The vertical axis points 70° counterclockwise from 
north. 

An eccentric disk may explain other features observed by HST (L93; King et al. 1995): 
(a) the region between PI and P2 is darker than would be expected from the superposition 
of two elliptical components (there may be a vacant region inside the nuclear disk); (b) the 
PI isophotes are twisted (isophote twists are a natural consequence of an eccentric disk 
whose eccentricity or thickness varies with radius); (c) at the highest surface brightness 
levels, both PI and P2 appear elongated along the P2-P1 axis, a natural feature of a disk 
viewed at a high inclination angle; (d) P2 exhibits a weak cusp and enhanced ultraviolet 
emission at its center, as might be expected close to a BH, while PI does not, as expected 
for a disk. 

The center of the bulge should coincide with the center of mass of the BH and the disk 
rather than precisely with the BH (i.e. P2). Thus in the best-fit model the center of the 
bulge should be displaced from P2 towards PI by 0”02. King et al. (1995) estimate that 
the center of the bulge is displaced from P2 towards PI by 0”05, while L93 estimate that 
the displacement is between 0 and 0”05. 

Figure 3 shows the mean velocity and velocity dispersion as a function of distance from 
P2"^. The left-hand panels show data from Kormendy (1988) along position angles PA= 55° 
(filled squares) and 38° (open squares), along with predictions obtained by convolving the 

The comparison of theoretical dispersion prohles with measured dispersions can be misleading, espe¬ 
cially when the light comes from stars with a range of velocity dispersions (the cool nuclear disk, the hot 
bulge, and the even hotter cusp of bulge stars near the BH). This complication is neglected in the present 
paper (see Kormendy 1988 and van der Marel 1994 for discussion). 
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3 . Rotation velocity and velocity dispersion as a function of distance from P2 in the best-fit 
model. The prohles have been “observed” with a Gaussian point-spread function of standard 
deviation cr*, through a slit of width s at position angle PA. The curves in the left-hand panels have 
parameters that approximate Kormendy’s (1988) observations: the solid curves have cr* = 0!^52, 
5 = o!^5, PA= 55°, and the dotted curves have cr* = 0!^65, s = l!^0, PA= 38°. The corresponding 
observations are shown as hlled and open squares, respectively. Note that positive r in Kormendy 
(1988) corresponds to negative x, i.e. PI is at x ~ — 0!^5. The prohles in the right-hand panels 
have ( 7 * = 0 !^ 37 , s = 0!^5, PA= 54°, corresponding to Bacon et al.’s (1994) observations which are 
shown as hlled circles. The dashed lines show the rotation and dispersion curves that would be 
observed with perfect resolution. 

best-fit model with the appropriate point-spread functions and slit widths. The profiles 
that would be seen at perfect resolution are shown as dashed lines. The right-hand panels 
show rotation and dispersion data from Bacon et ah (1994). 

First consider the comparison to Kormendy’s data. The model correctly predicts the 
approximate shape of the rotation profile, and the location of the zero-rotation point close 
to P2. The maximum rotation speed is too low by ~20%, suggesting that the model 
underestimates the prominence of the nuclear disk. The model correctly predicts the 
shape of the dispersion profile—although the bulge mass-to-light ratio and the mass of 
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the central BH have been chosen to optimize this fit. There are slight asymmetries in the 
model rotation curve: the maximum of the rotation curve on the PI side {x < 0) is at a 
larger distance than on the anti-Pl side ( — 1” compared to +0”8), and the zero-rotation 
point is displaced 0”1 towards PI from P2. There is no clear asymmetric signal in the 
data, although asymmetries would be masked by Kormendy’s decision to define the origin 
as the zero-rotation point of the rotation curve. The dispersion profile of the model is 
asymmetric, with dispersions on the anti-Pl side {x > 0) being larger than dispersions 
on the PI side; this agrees with the sense of the asymmetry in the data although the 
asymmetry in the model is smaller. 

The right-hand panels compare the model to data from Bacon et ah (1994). The model 
correctly predicts the amplitude of the rotation curve; the principal discrepancies are that 
(a) the linear part of the model rotation curve is too steep; (b) the zero-rotation point of 
the model curve is displaced from P2 towards PI by about 0”2, while the analogous point 
in the data is displaced away from PI by 0”l-0”2. The model dispersion profile correctly 
predicts that the dispersion is enhanced on the anti-Pl side although the width of the 
dispersion peak is too small; as in the case of the rotation curve, the agreement would be 
improved if the model curve were shifted away from PI. These discrepancies are probably 
minor given the spatial resolution of the data (sampling of 0”39 and seeing FWHM of 
0"87). 

All of these kinematic comparisons are with a single model that was designed mainly 
to match the photometry along the P2-P1 axis. No kinematic information was used in the 
fit except that the mass of the central BH and the bulge mass-to-light ratio were adjusted 
to match the dispersion curve. Simultaneous fits of eccentric disk models to both the 
photometry and the kinematics have not been explored. 


3 Discussion 

All eccentric disk models that match the sky brightness distribution should share several 
features: 

(a) If the disk lies in the same plane as the M31 disk (inclination 77°) then it cannot be 
too thin, or else the isophotes surrounding PI would be strongly flattened towards the 
P2-P1 axis. To fit the observations, the rms disk thickness h must exceed roughly 0.3r 
where r is the radius. In fact a thick disk is also required by dynamical arguments. 
The thickness grows on the two-body relaxation time, which is given by 
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where P = 27r(r^is the orbital period, Ar is the radial width of the disk, m 
is the typical stellar mass, Md is the disk mass, and In A is the Coulomb logarithm. 



Inserting approximate values for the M31 disk, 
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Thus the disk thickness h must exceed 0.3r in order that the relaxation time is longer 
than the Hubble time. 

(b) The surface brightness of the disk may be written /(r, cf)) where r is the radius in the 
disk and (f) is the azimuthal angle. Suppose that the disk is cold, the streamline with 
semi-major axis a has eccentricity e(a) and periapsis angle zu{a), the luminosity in 
the interval [a, a -\- da] is £(a)da, and the azimuthal origin is chosen so that zu = 0 
or TT. Then e and ru may be replaced by the single variable k(a) = e(a) cos w(a) 
( — 1 < k < 1), and 


T/ '^(®) 1 — k 

I(t 0 ) =- 

’ 27r(l — 1 — k — adk/dad 

. . £{a) _ l + k 

’ 27r(l — k'^y!'^ 1 + k -\- oAk/dad 
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where a is determined implicitly from r{a, /i = 0) = a[l — k{a)] and r{a, /i = vr) = 
a[l -|- k{a)]. The condition that streamlines do not intersect is \k -\- adkjda\ < 1. The 
surface brightnesses at a given semi-major axis on opposite sides are in the ratio 

.. adk / da 

/[r(a, 7r),7r] ^ ^ ~ 1-k .gx 

/[r(a, 0),0l 1 , adk/da' 

l + k 


Now let /i = 0 correspond to the direction from P2 to PI. Since the surface brightness 
at PI is higher than anywhere on the anti-Pl side, the right side of equation (6) must 
be less than unity, which requires dk/da > 0. Thus, for example, if PI coincides 
with the apoapsides of the disk streamlines, then ru = vr. A; < 0 and the eccentricity 
gradient dejda must be negative in the brightest parts of the disk (consistent with 
the monotonic decrease of e with increasing a in the ringlets of equation 2). 

(c) The density distribution along a Keplerian orbit can be expanded as a Fourier se¬ 
ries in azimuth. If the eccentricity is small, the dominant azimuthal wavenumber is 
m = 1. The asymmetry associated with PI has substantial components at azimuthal 
wavenumbers greater than 1; therefore the eccentricity of the disk in the region con¬ 
tributing to PI must be substantial, at least 0.3 and probably greater. 

A mechanism is needed to maintain apsidal alignment against differential (retrograde) 
precession caused by the bulge (in the best-fit model, differential precession between near¬ 
circular orbits at 0”5 and 1" is ~ 1 radian/Myr). One possibility is that the apsides are 
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trapped at a preferred orientation by the triaxial potential from the bulge (Stark 1977). 
However, for modest levels of triaxiality only highly eccentric orbits are trapped; also, if 
the disk is composed of stars on trapped orbits with periapsis angle ru, then we must 
still explain why the equivalent trapped orbits with periapsis angle zu -\- tt are not also 
populated. 

A more promising possibility is that the alignment is maintained by the self-gravity of 
the disk. The eccentric distortion can be regarded as a discrete nonlinear eigenmode of a 
self-gravitating disk orbiting in the external field of the BH and the bulge. The eigenmode 
has a pattern speed 0^, the common apsidal precession rate, which is much less than the 
typical Keplerian angular speed O = (GM,because the disk mass is much less than 
the BH mass. The properties of such disks have not been investigated much^, and we would 
like to know under what circumstances they exist, and what are their shapes and pattern 
speeds. 

A related problem is why the disk should be eccentric. Gurzadyan & Ozernoy (1979) 
and Syer & Clarke (1992) have suggested that eccentric accretion disks form through tidal 
disruption of an object passing by the BH (most likely a giant molecular cloud). Another 
possibility is that the eccentricity arises from an unstable m = 1 normal mode of the 
axisymmetric disk. A third mechanism, perhaps the most promising, is the influence of 
dynamical friction from the bulge on a nearly circular disk. Suppose for the moment that 
the bulge does not rotate and its velocity dispersion tensor is isotropic, so that dynamical 
friction removes energy from the disk, E < Q (Tremaine & Weinberg 1984). Dynamical 
friction also removes angular momentum, at a rate J = EjQp. Thus if 0 < Op ^ O, 
angular momentum is removed much faster than energy so the disk orbits become more 
eccentric: in other words there is a secular instability that causes the eccentricity of the 
disk to grow. On the other hand, if Op < 0 the friction adds angular momentum so 
the disk becomes less eccentric. Disks can have positive or negative pattern speeds: the 
bulge potential causes retrograde precession but the self-gravity of the disk causes prograde 
precession. Since the mass of the disk is larger than the mass of the bulge inside the disk, 
most or all eigenmodes have Op > 0 so their amplitudes grow. If the bulge rotates with 
mean angular speed O;,, then only eigenmodes with Op > O;, will grow. An unsolved 
problem is what limits the growth of the secular instability. 

The presence of an eccentric disk is difficult to reconcile with the arrival of 10®'^ Mq 
BHs every ~ 10® yr, as required if the halo is composed of such objects (Lacey & Ostriker 
1985; Xu & Ostriker 1994), since each inspiraling BH would destroy the disk. 

The model described here implies that the size of the dark object at the center of M31 
should be < 0”25 ~ 1 pc (the periapsis distance of ringlet 1) so that the force held seen 
by the eccentric disk is nearly Keplerian—otherwise the relatively weak self-gravity of the 
disk could not maintain apsidal alignment. This upper limit is more stringent than limits 
based on variation of mass-to-light ratio with radius (Richstone, Bower & Dressier 1990) 
and close to the lower limit for the size of a stellar cluster, which must have half-mass 
radius > 0”1 (Goodman & Lee 1989). 


^ Except for the case of narrow cold disks with sharp edges (Goldreich & Tremaine 1979). 
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Higher resolution spectroscopic observations, improved photometry from HST, and 
measurements of the shape of the line-of-sight velocity distribution can all provide sensitive 
probes of eccentric disk models in the near future. Equally important tasks are to develop 
a better theoretical understanding of the dynamics of collisionless eccentric disks, and 
self-consistent models of eccentric disks that predict the line-of-sight velocity distribution. 

I thank my HST collaborators Ed Ajhar, Yong-Ik Byun, Alan Dressier, Sandra Eaber, 
Karl Gebhardt, Carl Grillmair, John Kormendy, Tod Lauer, and Doug Richstone for many 
conversations and insights about galactic cores and nuclei. I also thank Roland Bacon, 
Gerry Quinlan, Jerry Sellwood, Jerry Ostriker, and Dave Syer for discussions about the 
nucleus of M31, and Jerry Ostriker for introducing me to this interesting system. This work 
was carried out during a visit to the Institute of Astronomy, Cambridge, which was made 
possible by assistance from the Killam Program of the Canada Council and the Raymond 
and Beverly Sackler Eoundation. The research was also supported by grants from NSERC. 
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